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ABSTRACT 



43 

O t ' The linewidth (a) - size (R) relationship of the interstellar medium (ISM) has 

I \ been extensively measured and analysed, in both the local ISM and in nearby normal 

galaxies. Generally, a power-law describes the relationship well with an index ranging 
from 0.2 - 0.6, and is now referred to as one of "Larson's Relationships." The nature 
of turbulence and star formation is considered to be intimately related to these rela- 
tionships, so evaluating the cr — R correlations in various environments is important for 
developing a comprehensive understanding of the ISM. We measure the linewidth-size 
relationship in the Central Molecular Zone (CMZ) of the Galactic Centre using spec- 
tral line observations of the high density tracers N 2 H+, HCN, H 13 CN, and HCO+. We 
construct dendrograms, which map the hierarchical nature of the position-position- 
velocity (PPV) data, and compute the linewidths and sizes of the dendrogram-defined 
structures. The dispersions range from ~ 2 - 30 km s _1 in structures spanning sizes 2 - 
40 pc, respectively. By performing Bayesian inference, we show that a power-law with 
exponent 0.3 - 1.1 can reasonably describe the cr — R trend. We demonstrate that the 
derived cr — R relationship is independent of the locations in the PPV dataset where a 
and R are measured. The uniformity in the cr — R relationship indicates that turbulence 
in the CMZ is driven on the large scales beyond ^ 30 pc. We compare the CMZ cr — R 
relationship to that measured in the Galactic molecular cloud Perseus. The exponents 
between the two systems are similar, suggestive of a connection between the turbulent 
properties within a cloud to its ambient medium. Yet, the velocity dispersion in the 
\ CMZ is systematically higher, resulting in a scaling coefficient that is approximately 

five times larger. The systematic enhancement of turbulent velocities may be due to 
the combined effects of increased star formation activity, larger densities, and higher 
pressures relative to the local ISM. 

Key words: ISM: clouds - ISM: molecules - ISM: structure - stars: formation - tur- 
bulence 
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1 INTRODUCTION 

Turbulence is observed on all scales larger than the size 
of the densest star forming cores (~0.1 pc), and is con- 
sidered a key process influenc ing star formation i n the 
interstellar medium (IS M, see iMac Low &; Klessen 2004; 
iMcKee &; Ostrikei] 120071 and references therein). Atomic 
and molecular line observations provide detailed informa- 
tion about the morphological, thermal, chemical, and dy- 
namical state of the ISM, such as the masses, velocities, 



and sizes of a variety of structures, from the densest star- 
forming cores, larger molecular clouds ( MCs), and the sur- 
rounding, volume filling diffuse gas (e.g. lYoi mg & Sc ovillel 
Il99ll : iKalberla &; Kerdl2009l : iFukui fe Kawamurall201(j and 
references therein). Accurately measuring these properties 
and the relationships between them is necessary for devel- 
oping a complete understanding of turbulence in the ISM, 
including its role in the star formation process. 

Spectral line observations of the ISM show strong 
power-law correlations between the mass M and linewidth a 
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with their projected size R. Observations of Milky Way MCs 
generally exhibit exponents ^2 and ~ 0.2-0.6 for th e mass- 



size and linewidth-size relationship, resp ectively (ILarson 



1981 


; Solomon et al. 1987; Leisawitz 1990; Falgarone et al. 


1992 


; Caselli & Mvers 1995; Kauffmann et al. 2010a, b; 


Beaumont et al.ll2012f). Similar trends are observed for ex- 



tragalactic clouds (jBolatto et al.1 12008 : ISheth et al.l l2Q08h . 
These results are interpreted as evidence that MCs are vit- 
alized objects, and are generally referred to as "Larson's 
Relationships." 



That the M — R and a — R relationships portray sim- 
ilar power law scalings across a range in environments sug- 
gests that the underlying p hysical processes drivi ng ISM dy- 
namics is "universal" (e.g. iHever &; Bruntll20o3 ). However, 
it is not clear if this "universality" applies to systems signifi- 
cantly different than the Milky Way or Local Group GMCs, 
as it has not been possible to sufficiently resolve the ISM 
in external systems. Measuring these relationships in vary- 
ing environments, such as (ultra) luminous infrared galaxies 
([UjLIRGs) or similar starbursting regions like the Galactic 
Centre, is required to comprehensively understand the ISM 
throughout the universe. 

To understand the nature of turbulence in a starburst- 
ing environment, we investigate the linewidth-size relation- 
ship of high density structures in the Central Molecular Zone 
(CMZ) of the Galactic Centre (GC). The CMZ is denser 
than the ISM outside the GC, with densi ties comparable to 
the s tar forming cores within MCs (see [Morris &; Serabvnl 
1996, and references therein). The star formation rate in the 
CMZ is ~1.5 orders of magnitude higher than the rate mea - 
sured in the solar neighborhood (|Yusef-Zadeh et all l20Q9h . 
Therefore, the star formation activity in the GC is in many 
ways similar to starbursting systems like LIRGs. Given its 
relatively close proximity (~ 8 kpc), the GC affords us the 
best opportunity to resolve the structure of the ISM in star- 
bursting environments, and probe the nature of the dynamic 
ISM. 

In this work, we measure the linewidth-size relationship 
in t he CMZ using Mopra obs ervat ions of dense ga s trac - 
ers. iMivazaki fe TsuboJ (|20Q0h and lOka et al.l (|l998l . l200lh 
present an analysis of the Larson scaling relationships from 
Nobeyama observations of CS and 12 CO, respectively. They 
show that the a — R relationship in the CMZ is similar 
to that in the Milky Way disk, but with enhanced veloci- 
ties. We extend their analysis by considering four high den- 
sity tracers, N2H + and H 13 CN with low opacity, as well 
as HCN and HCO + with high opacity. Further, the Mopra 
coverage of the CMZ is fully (Nyquist) sampled, allowing 
for thorough and systematic identification of dense clouds. 
By constructing hi erarchical structure trees, or dendrograms 
(jRosolowskv et al.ll2008h . we identify contiguous features in 
the observed position-position-velocity (PPV) cubes of the 
CMZ. One key difference between our study and previous 
efforts is the manner in which "clouds" are identified. Den- 
drograms map the full hierarchy of structures, such that 
ISM structures which are not necessarily part of the densest 
clouds are included in our analysis. We assess the linewidth- 
size relationship of the ensemble of these structures, and 
compare it to the well known a — R trends found in the 
local ISM. 

This paper is organised as follows. The next section de- 
scribes the observations of high density tracers in the CMZ. 



In Section E3J we describe the method we use to iden- 
tify contiguous structures in the observed data. Section [3] 
presents the linewidth-size relationship. In Section [H we dis- 
cuss how the measured a — R relationship relates to the un- 
derlying turbulent velocity field, and compare it to the trend 
in the Milky Way molecular cloud Perseus. We summarise 
our findings in Section [5] 



2 OBSERVATIONS AND ANALYSIS 
TECHNIQUE 

2.1 Mopra survey of the CMZ 

In our analysis of the linewidth-size relationship of the CMZ, 
we use the observations from the 22 metre Mopra millimetre- 
wave telescope, which is operated by the Australia Telescope 
National Facility (ATNF). The Mopra CMZ survey covers 
the 85-93 GHz window, covering 20 spectral lines including a 
number of dense gas tracers like the HCN, HNC and HCO + 
mole cules, as well as the cold core species N2H + (|jones et al.l 
The total area covered by the survey is 2.5x0.5 deg , 
with 40" spatial resolution (~2 pc at a distance of 8 kpc), 
and spectral resolution < 0.9 km s _1 . For our analysis, we 
use smoothed datasets with 3.6 km s _1 spectral resolution 
of observations of the lowest rotational transitions of four 
molecules: N 2 H+ (rest frequency 93.173 GHz), HCN (88.632 
GHz), H 13 CN (86.340 GHz), and HCO+ (89.189 GHz). Fur- 
ther information about the Mopra CMZ survey, including 
reduction tech niques, inten s ity m aps, and velocity structure, 
is described in I Jones et al 

2.2 13 CO FCRAO observations of the Perseus 
molecular cloud 

In Section [4.41 we compare the linewidth-size relationships 
in the CMZ to the Milky Way molecular cloud Perseus. We 
use the 13 CO J= l-0 line (110 201 G Hz) data from the COM- 
PLETE Survey (jRidge et al.ll2006h . The data are approxi- 
mately Nyquist sampled, with a pixel size of 23" (.028 pc at 
a distanc e of 250 pc) , and h ave a velocity resolution of 0.066 
km s- MRidge et al.1 (|2006h describes the COMPLETE ob- 
servations of Perseus, including 12 CO and 13 CO, as well as 
HI, dust extinction and emission maps of the cloud. 

2.3 Measuring the linewidth-size relationship 

In order to measure the linewidth and sizes of structures, 
we first need to identify the relevant structures using a well 
defined criterion. One rather simple definition of a struc- 
ture is any contiguous region with (approximately) constant 
density. Unfortunately, we only have information about the 
intensity of a number of spectral lines at a range of velocities, 
along lines of sight (LoS) through the CMZ. As the CMZ 
has extremely dense structure throughout its volume, veloc- 
ity integrated intensity maps would merge a number of dense 
features along the LoS, thereby smoothing out the intrinsic 
struc ture of the ISM i n the resulting 2-dimensional (2D) map 
(e.g. Pichardo et al.l |2000| ; iBallesteros-Paredes &; Mac Low! 

1 All data from the Mopra CMZ survey is publicly available at 
http:/ /www. phys.unsw.edu.au/mopracmz 
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Figure 1. Dendrogram of a hierarchical cloud. Left: simple cloud 
with three levels of hierarchical substructure. Right: Dendrogram, 
or structure tree, of the cloud. Each distinct (closed contour) re- 
gion, is identified as a structure (vertical lines). Structures on the 
top of the tree are "leaves." The structures which enclose multi- 
ple higher level contours are "branches," and the structure which 
encloses all other structures is a "root." 



120021 : iGammie et alJl2003l : IShettv et al.ll2Q10h . Accordingly, 
contiguous regions in the 3-dimensional (3D) PPV space 
should be able to more accurately identify distinct struc- 
tures. As the ISM contains hierarchical structure, direct 
clump decomposition methods may not accurately identify 
different levels of substructure, especially in high density, 
turbulent environments. We thus employ the de ndrogram 
decomposition technique ([Rosolowskv et all [2008) to iden- 
tify and characterize contiguous structures from the PPV 
cubes of the CMZ. 

We provide a brief overview of dendrog rams here, and 
refer the reader to iRosolowskv et all (|2008h for a more de- 
tailed description of the structure decomposition algorithm. 
Dendrograms provide a "structure tree" which quantifies the 
hierarchical nature of any dataset. Consider the cloud de- 
picted in Figure [1] The cloud is composed of three different 
levels, or intensities, with the base level (light blue struc- 
ture) enclosing all higher levels. Each structure at a given 
level is distinguishable from adjacent structures if a closed 
iso-intensity contour can enclose the whole structure. As the 
intensities of the features do not merge in a contiguous man- 
ner, they are identified as separate structures. Both struc- 
tures in the second level encloses higher intensity regions, 
which are themselves separated into two distinct structures. 

The structure tree in Figure [T] maps the hierarchical na- 
ture of the model cloud. Each vertical line corresponds to a 
distinct structure in the cloud. Those on the top of the tree 
are referred to as "leaves." The lowest level is the "root" 
of the tree. Any level which encloses multiple structures re- 
quires "branches" to accommodate the hierarchy. The leaves 
at top of the tree shows the highest density structures which 
do not enclose any further substructure. 

The mapping illustrated in Figure [1] can also be per- 
formed for 3D data, such as a PPV cube. Each structure 
corresponds to an iso-intensity region in PPV space. Ad- 
ditional information can be derived from the dendrogram- 
identified regions. As the structures from a PPV cube de- 
marcate iso-intensity regions, the velocity dispersions a of 
this structure may be directly computed. We compute the 
intensity weighted velocity dispersion of each structure 
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Figure 2. Linewidth-size relationship in the CMZ, as measured 
within dendrogram identified structures in N2H+, HCN, H 13 CN, 
and HCO+. Filled symbols correspond to "leaves" which do not 
enclose additional higher level structures. Open symbols are struc- 
tures which do contain higher level structures. Lines show the best 
(X 2 ) power-law fits. 
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where the summations are taken over all zones of the iden- 
tified structures, I v is the line intensity, v is the LoS ve- 
locity, and and v is the intensity weighted mean velocity, 
Yl(Iyv)/y^ Iv We also define the size R of each struc- 
ture using the projected area, R = (N pro jAxAy/7v)2, where 
iVproj is the number of projected pixels, and AxAy is the 
area of the resolution element. 

To ensure that the dendrogram denned structures are 
well resolved, we exclude any features with R < 2 pc, and 
a < 2 km s _1 . Since we are interested in the dispersion of the 
contiguous and distinct structures, and not the whole ISM 
in the CMZ, we also only consider features with R < 40 pc. 
Such a size limit may further minimise contamination from 
the superposition of features along the LoS. Finally, we also 
exclude any region which occurs within three voxels (zones 
in a PPV dataset) from the boundary of the PPV cubes. 
Under these conditions, we minimise the likelihood of iden- 
tifying features dominated by observational noise, resulting 
in well denned and contiguous structures PPV space. We 
have adopted a uniform "dendrogramming" strategy such 
that each structure contains at least 50 voxels. This con- 
dition ensures that the dendrogram identifies real extended 
structures in the PPV cube, which are unlikely to be purely 
noisy features. From the measured value of a and R, we can 
now assess the linewidth-size relationship of the ensemble of 
dendrogram denned structures. 
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3 THE LINEWIDTH-SIZE RELATIONSHIP IN 
THE CMZ 

Figure [2] shows the a — R relationship of the dendrogram de- 
fined structures from the N 2 H+, HCN, H 13 CN, and HCO+ 
observations of the CMZ. The largest structures identified 
in each tracer all have dispersions of 20-30 km s _1 , which 
is equivalent to the global linewidth of the CMZ. These 
structures generally enclose higher level, brighter features. 
Towards smaller scales, the dispersions systematically de- 
crease, reaching ^2 km s _1 in structures with R~2 pc. The 
open symbols in Figure [2] are regions which contain at least 
one higher level structure. The closed symbols correspond 
to the leaves on the top of the structure tree (see Fig [T]) . 
The ensemble of structures exhibit a systematic decrease in 
linewidths from the largest scale structures to the smallest, 
brightest features in all the observed lines. 

That the structures in all tracers span a limited region 
in a — R space, as well as depict similar trends, provides con- 
fidence that the apparent a — R relationship in Figure [2] is 
real. N2H + has a critical density of ^10 5 cm -3 , and is a good 
tracer of cold, dense cores. HCN and HCO + tra ces more ex- 
tend ed dense gas in the ISM (see Figure 2 in Ijones et al.l 
2012), with similar critical densities for each J=l-0 emis- 
sion lines. They can be moderately op tically thick in the 
CMZ, as analysed bv I Jones et al.l (|2012l ). Hence, we perform 
the analysis on the weaker, but optically thin H 13 CN line, 
obtaining essentially the same results; this indicates that 
optical depth has not significantly biased the trends. Both 
HCN and N2H + exhibit hyperfine splitting, being split over 
^14 km s _1 into 3 groups of lines in the ratio 3:5:1. Though 
this is less than the typical turbulent width for the gas in 
the CMZ, multiple smaller scale dendrogram structures in 
N2H + and HCN may actually be part of a single feature. 
On the other hand HCO + , which is not split into any sub- 
components, produces very similar results, indicating that 
the derived linewidth-size relation of the ensemble is not 
significantly affected by line splitting. 

There are some differences between the linewidth-size 
trends from the various tracers shown in Figure El Most 
notably, the HCN and HCO + have a number of larger scale 
structures with R ^ 3 pc and a £ 5 km s _1 that do not 
occur in N2H + . These differences can be largely attributed 
to the relative extents of the different tracers, and for some 
structures the foreground absorption of HCN and HCO + 
emission. 

For some of the large R and low a points, a closer inspec- 
tion of these structures reveals that they are further divided 
into separate features in the N2H + dataset. Figure [3] shows 
an example of one such structure in N2H + and HCN. The 
HCN dendrogram-defined contour corresponds to the leaf 
with the lowest a— 2.1 km s _1 , with R— 5.1 pc. The dendro- 
gram of the N2H + feature, on the other hand, is more com- 
pact with R— 2.8 pc and a— 2.8 km s _1 . N2H + traces cold, 
dense cores, and as it can be destroyed by reactions with 
CO, it is expected to be present in the den sest and most qui- 
escen t regions where CO is depleted (see IBergin &; Tafallal 
2007, and references therein). N2 H + is also found to be les s 
extended than HCN and HCO + (jPurcell et al .11 20061 . liocih . 
In addition, the HCN and HCO + line emission has higher 
optical depth. Hence the dendrogram of an HCN and HCO + 
structure may therefore demarcate a larger region surround- 



ing a dense peak than does N2H + , as shown in Figure [3] 
Consequently, an extended HCN structure may be found to 
break up into a number of compact N2H + features, with 
similar measured linewidths. The discrepancies in the a and 
R provided by the different tracers may not necessarily re- 
flect morphological variations of the gas, but rather diversity 
in chemical abundances. The subdivision of features in the 
N2H + observations results in a more uniform a — R relation- 
ship, and correspondingly a better defined trend for a larger 
range in sizes. 

As discussed by Ijones et al.l (|2012l ). there are strong 
foreground absorption features in the HCN and HCO + 
datasets, possibly due to the presence of the spiral arms 
along the LoS. The sharp drop in intensity of some of the 
large scale features across the velocity channels correspond- 
ing to -52, -28, and -3 km s _1 leads the dendrogram al- 
gorithm to identify multiple structures. Consequently, such 
an extended feature only spans a limited range in velocity, 
thereby resulting in a structure with large R but low a. The 
division of large scale structures due to foreground absorp- 
tion along the LoS does not occur as much in N2H + and 
H 13 CN, as those molecules are not as pervasive, resulting in 
substructure that can be more easily differentiated from its 
background, as described above. 

The a — R trends shown in Figure [2] indicate that there 
is a large range of possible dispersions for a given size, 
reaching ^ 7 km s _1 in the HCN dataset. Several fac- 
tors can affect the amount of substructure and measured 
properties from any observation, including signal-to-noise 
ratio, resolution, chemical abundance inhomogeneities, and 
radiative transfer effects (opacity, sensitivity to tempera- 
ture and density). Further, local systematic effects may in- 
fluence some regions of the ISM more than o thers. For in- 
stance, |llodrigu£z^^r^^ (|2008l ) suggest that 
infalling gas along spiral arms and through the large cloud 
complex at longitude 1.3° lead to the asymmetry of the 
CMZ. We find that the measured linewidths in this region 
are not particularly discrepant from the ensemble, as the 
values of a and R measured in that complex fall within the 
range provided by the other structures. The intrinsic scat- 
ter of the measured a for structures with similar sizes may 
be reflective of the differences in environmental conditions 
across the ISM. 

Despite the numerous physical processes that likely 
influence the gas dynamics in the CMZ, Figure [2] shows 
that the linewidths are correlated with the sizes. The Pear- 
son's correlation coefficient of the N 2 H + , HCN, H 13 CN, and 
HCO + linwidth-size data (in log space) are 0.84, 0.58, 0.86, 
and 0.72, respectively. The ensemble may thereby be de- 
scribed as power-law relationships. To quantify the power- 
law, we fit 

in log space to the dendrogram derived linewidths and sizes. 
The fit estimates the power law index b and coefficient A. 
We first employ a x 2 fit, and in order to comprehensively 
account for the uncertainties, we also perform a Bayesian 
linear regression analysis. 
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Figure 3. A single a) N2H + and b) HCN channel map (corresponding to velocity Vlsr ~140 km s _1 ) showing a dense filamentary 
region near to the Galactic Centre. White contours show the cross-section of one dendrogram defined iso-intensity surface at this velocity. 
The linewidths and sizes of the whole dendrogram structures are 2.8 km s _1 and 2.8 pc in N2H+, and 2.1 km s _1 and 5.1 pc in HCN, 
respectively. 



Table 1. Best (x 2 ) fit parameters of linewidth-size relationship 



Tracer 


N 2 H+ 


HCN 


H 13 CN 


HCO+ 


Power law index b 


0.67 


0.46 


0.78 


0.64 


Coefficient A 


2.6 


3.8 


2.6 


2.1 



1 The formal la errors in b and A are all ^ 0.06 and 1.2, 
respectively. 

3.1 Least squares (x 2 ) fit 

As a first estimate of the power- law index and coefficient, we 
apply a standard \ 2 fit- Each panel in Figure [2] shows the \ 2 
fit lines, and Table [T] provides the best fit parameters. The 
maximum error in the best fit slopes is 0.06, obtained from 
the HCO + structures. This error is the formal x 2 error, and 
in this case is primarily driven by the large number of points 
which results in rather low values. The resulting error esti- 
mates suggest that to a high degree of probability (3a) the 
N2H + and HCN features have different slopes. In principle, 
it is possible that different molecules exhibit different a — R 
slopes, since they may span different regions in the ISM. 
Such disparity may be expected when considering tracers of 
different densities (e.g. 12 CO vs. HCN), which would cor- 
respond to regions with variable extents. However, in this 
study we are focusing only on high density tracers, and a 
direct comparison of the structures in the PPV datacubes 
do not show obvious signs that the tracers span regions of 
variable spatial extents. Additionally, Figure [2] suggests sig- 
nificant overlap of a and R for each of the tracers. Therefore, 
the low error estimate in the % 2 fit may not accurately char- 
acterize the uncertainty in the linewidth-size relationship. In 
order to rigorously account for the errors, as well as ascer- 
tain the full range of possible best fit parameters, we employ 
a Monte Carlo analysis and a Bayesian linear regression fit 
for estimating the coefficient and power-law index. 



3.2 Uncertainties in the dendrograms 

Dendrograms identify contiguous structures in the PPV 
cubes based on the intensity at each position and LoS ve- 
locity. Any uncertainties in the intensities may lead to inac- 
curate identification of structures, and thereby an erroneous 
a — R relationship. To assess the magnitude of these uncer- 
tainties, we perform a simple bootstrapping analysis. 

For this analysis we use the N2H + dataset, which has a 
well denned a—R relationship. We estimate a la uncertainty 
0.04 K from a region far from any significant emission. For 
each realization in the Monte Carlo simulation, we modify 
each voxel in the PPV cube by an additive factor e, where e 
is drawn from a normal distribution with mean and 0.04 
K standard deviation. After constructing the structure tree, 
we fit the linewidths and sizes as described above. 

After performing this exercise 100 times, we measure 
a standard deviation in the best fit parameters b and A to 
be ^ 0.01. This value is lower than the errors returned 
by the \ 2 fit of the original unperturbed dataset. This test 
indicates that the variance in the number of structures, as 
well as the uncertainty in the estimated linewidths and sizes 
due to noise, hardly influences the fit a — R relationship. 
We have performed a similar analysis on the other three 
datasets, and as with the N2H + analysis, we find that the 
range in the best fit parameters to be smaller than the errors 
computed in the \ 2 fit- We can thus be confident that the 
dendrograms of the original dataset are not sensitive to noise 
in the PPV datacubes. We now turn our attention to the 
affect of resolution, for which we employ a Bayesian linear 
regression analysis. 



3.3 Bayesian regression fit 

One of the primary advantages of Bayesian inference is that 
the uncertainties in the measur ed data are rigorously and 
self-c onsistently treated (e.g. Gelman et al. 2004; Krusc hke 
l201lh . A Bayesian fit considers the error in each measured 
quantity, in our case a and R, to be drawn from some 
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a priori defined distribution. The choice of the distribu- 
tion should reflect the uncertainties in the measurement. 
Commonly, Bayesian inference utilizes Markov Chain Monte 
Carlo (MCMC) routines to sample the probability distribu- 
tion of the fit parameters, given the measured data and the 
assumed uncertainties. The result of the Bayesian inference 
is a joint probability — or posterior — distribution of the re- 
gression parameters. This probability distribution function 
(PDF) accounts for the assumed measurement error distri- 
butions, and therefore provides well defined uncertainty es- 
timates on each of the fit parameters, as well as the corre- 
lations between them. We refer the reader to iKellvl ([2007) 
for details of the Bayesian fitting algorithm we use in our 
analysis 

For our estimate on the uncertainties in the size we use 
the effective spatial resolution of the observational datasets. 
For the error estimate in the velocity dispersion, we choose 
a conservative value of 2 km s _1 , which is about twice the 
resolution of the original observations. Note that though the 
resolution of the smoothed datacube (3.6 km s _1 ) is larger, 
any given contiguous region in a PPV cube may produce 
a dispersion lower than the effective resolution, as the dis- 
persion is simply a measure of the deviation from the mean 
value. The smoothing primarily affects the very smallest fea- 
tures with narrow linefl which are excluded from our anal- 
ysis given our criteria described in Section 2.3. 

In deriving the linewidths and sizes, we extract contigu- 
ous regions in PPV space defined by the dendrogram algo- 
rithm. Though this process introduces additional sources of 
error, we have shown in Section 13.21 that these uncertain- 
ties are minimal, and hardly affect the a — R relationship. 
Therefore, in our Bayesian fit we only consider the uncer- 
tainties due to resolution (e.g. beam smearing). We are thus 
effectively assuming that resolution is the dominant source 
of uncertainty, and that the errors in log(cr) and \og(R) are 
normally distributed, with dispersions 2 km s _1 and 2 pc, 
respectively. 

Figure 3] shows the Bayesian fit of the N2H + a — R 
relationship. The peak of the posterior, corresponding to 
a = 2.2i? ' 78 , is shown as the blue line. As expected, this 
line is similar to the \ 2 fit. The gray lines in Figure U are 
five random draws from the Bayesian posterior. The inset 
plot shows the marginal probability distribution of b. The 
95% (2a) highest density interval, or HDI, is marked by the 
vertical dashed lines, b £ 0.41 —1.13. This HDI provides a 
quantitative prediction of the range in possible b, given the 
uncertainties in the measurements of a and R. The distri- 
bution of b inferred by the Bayesian fit is six times larger 
than the uncertainty in the % 2 estimate. 

Figure [5] provides the marginal PDFs of the fit indices 
for all four tracers. The H 13 CN PDF is noticeably wider, re- 
flective of the lower number of identified structures in that 
tracer (Fig. [2]). Table [2] indicates the parameters from the 
peak of the posterior and in the 95% HDI for all the tracers. 
There is substantial overlap between the PDFs from the dif- 
ferent tracers, in both the coefficient and the exponent. It is 



2 IDL routines for the Bayesian linear regression fitting algorithm 
is publicly available at http://idlastro.gsfc.nasa.gov 

3 Line profiles along individual LoS are generally broader than 
3.6 km s _1 



N 3 H + 




1 10 
R (pc) 



Figure 4. Bayesian regression of linewidth-size relationship 
in N2H+ structures. Open circles are the dendrogram derived 
linewidths and sizes. The blue line corresponds to the fit with the 
highest probability. The gray lines show five random draws from 
the Bayesian posterior. The inset shows the marginal probability 
distribution P(b) of the power-law index. Vertical dashed lines in- 
dicate the 95% highest density interval (HDI), which ranges from 
0.41 - 1.13. 




2.0 



slope b 

Figure 5. Marginal distributions of the index b from the Bayesian 
regression analysis of the a — R relationship from N2H+ (thick 
black), HCN (thin red), H 13 CN (dashed green), and HCO+ (dot- 
ted blue). Table [2] provides the mode and the 95% highest density 
intervals of these PDFs. 



common practice in Bayesian inference to consider the full 
95% HDI as the reasonable range of fit parameters, rather 
than any single "best-fit" point estimate, so that the mea- 
surement uncertainties are taken into account. Under this 
framework, we cannot conclude that the evidence suggests 
a different a — R relationship for the various tracers. 
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Table 2. Results from Bayesian inference of linewidth-size parameters 



Tracer 


N 2 H+ 


HCN 


H 13 CN 


HCO+ 


Highest probability for b 


0.78 


0.62 


0.66 


0.79 


95% HDI of b 


0.41-1.13 


0.27-1.02 


0.21-1.66 


0.45-1.30 


Highest probability for A 


2.2 


2.7 


2.9 


1.7 


95% HDI of A 


0.9-3.8 


1.1-4.7 


0.0-5.8 


0.5-2.7 



4 DISCUSSION 

4.1 Do the growth of structures affect the a 
relationship? 



R 



The uniform a — R trend over the full range in spatial scale 
shown in Figure [2] suggests some connection between the 
densest structures and the larger scale ISM in the GC. As 
we discuss further in Section [4~4l this trend is also observed 
for clouds in the Milky Way disk, and provides evidence for 
large-scale turbulent driving. One closely related question is 
whether the growth of dense structures within a turbulent 
medium influences the underlying velocity field. 

To investigate the impact of dense clouds on the 
linewidth-size relationship, we perform a similar analysis to 
that described in Section [231 but with a few modifications. 
First, we assess the role of intensity weighting, simply by 
excluding the weighting factor I v in Equation [T] Next, we 
also measure the linewidth in random locations in the PPV 
cubes, using the dendrogram defined contours placed at dif- 
ferent locations throughout the PPV cubes. For the latter 
exercise, we shift the contours to a different random location, 
such that the contours no longer demarcate iso- intensity sur- 
faces. As a result, the shape of the contour is preserved, but 
the PPV region within the shifted contour may consist of a 
number of dense complexes, diffuse gas, or some combina- 
tion of both. To ensure that results are statistically robust, 
we perform numerous realizations of this exercise. Through 
these modified analyses, we can investigate the extent which 
the growth of dense structures and their internal processes 
influences the a — R relationship. 

Figure [6] shows the results from this modified analysis. 
The open circles are the original intensity weighted N2H + 
a — R points, and the solid line shows the fit derived from 
the maximum of the Bayesian posterior discussed previously 
and presented in Figured The smaller solid circles are the 
measurements where there is no intensity weighting. There 
is almost no difference between these two sets of points, since 
the intensity range within a given structure is limited. 

The red open triangles in Figure [6] show the results of 
one realization where the original dendrogram contours have 
been randomly shifted to a different PPV location. For this 
particular realization, the contours have been shifted by a 
minimum of 50 voxels in each of the three dimensions of 
the PPV cube, plus an additional random value drawn from 
a uniform distribution between and 20. There is still a 
clear trend of increasing linewidth with increasing size, with 
slightly more scatter than the original a — R distribution. 
The solid red line shows a = 2.1R 79 , which corresponds to 
the maximum of the Bayesian posterior for this dataset. This 
result is very similar to the original N2H + relationship, and 
the HDIs for the two PDFs almost completely overlap. We 
have performed ten different realizations of such a random 




Figure 6. Linewidth-size relationship of the N2H+ points. Black 
open circles are the original N2H+ measurements, and the black 
line is the most probable Bayesian inferred fit a = 2.2R 078 . Black 
solid circles show the linewidths computed without any intensity 
weighting. Red triangles are the measurements from random loca- 
tions in the PPV cube, using dendrogram defined contours. The 
red line is the most probable Bayesian inferred fit, a = 2.1R - 79 . 



shifting of the dendrogram contours, varying the amount by 
which the contours are shifted, and each one produces very 
similar results which are statistically indistinguishable from 
the original N2H + relationship. We have also performed this 
analysis on the HCN, HCO + , and H 13 CN datasets, and find 
the resulting a — R relationships to be statistically indistin- 
guishable from the corresponding original results. 

As clouds begin to grow, self-gravity and other internal 
process will begin to play a role in their subsequent dynam- 
ics. If internal processes dominate the ensuing dynamics of 
the clouds, we should find a distinct difference between the 
a — R computed from the dense structures, and that mea- 
sured from random regions irrespective of the presence of 
dense structures. That the measured a — R relationship can 
be reproduced regardless of the precise locations (in PPV 
space) where the linewidths are measured indicates that the 
dense cloud dynamics is predominantly governed by the un- 
derlying turbulent flow, with little contribution from internal 
sources of turbulence. 

In fa c t, a very similar analysis was carried out by 
>a et al. (199 (1), wh o concluded that the interpretation by 
Solomon et all (Il98<t ) that GMCs are bound arises due to 



the spatial and spectral crowding of clouds, in combination 
with the method employed to identi fy clouds in PPV cubes. 
Our analysis is in agreement with llssa et al ] (jl99d ). The 
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measured a — R relationship cannot fully reveal the dynam- 
ical state of the dense clouds. These results indicate that tur- 
bulence within dense clouds are driven on the largest scales 
in the ISM, a topic we further discuss in Section [4~4l 

Besides revealing that turbulence must be driven on the 
largest scales, the measured a — R relationships may delin- 
eate a lower limit to the turbulent driving scale. Applying 
different statistical tools on numeri cal simulations, such as 
struct ure functions and A- variance, Ossenkopf &; Mac Lowl 
(2002) showed that beyond the driving scale a will remain 
nearly constant with increasing R. Moreover, dendrograms 
are capable of identifying the flattening of the linewidth- 
size relationship towa rds the scale s at w hich turbulence is 
i nject ed, as shown by IShettv et al.l (|2010h and I Shettv et al.l 
(|201lL see their Fig. 13) using numerical simulations. Al- 
though the N2H + , HCN, and HCO + a — R pairs show some 
semblance of a decrease in b at the largest scales (Fig. [2}, 
there is no strong evidence for such flattening. Since a in- 
creases with R up to the largest scales in our analysis, 
we conclude that turbulence in the CMZ is predominantly 
driven on scales £ 30 pc. Note that though turbulence is 
driven on the large scales, it produces a velocity field which 
must be coherent on much smaller scales Quantifying this 
coherence may further reveal additional properties of the 
turbulent ISM. 
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Figure 7. Linewidth-size relationship of spherical (blue circles) 
and cubic (red squares) regions measured in the N2H+ PPV 
dataset. For comparison, the horizontal green dashes are the pre- 
dicted mean values of a from the Bayesian posterior (see Fig. HJ. 
The green vertical lines show the 95% HDI of the predicted values 
of a for each R. 



4.2 The relationship between the shape of dense 
structures and the a — R relationship 

So far, we have defined the contiguous structures in PPV 
space via dendrograms. That is, their shapes are governed 
by the morphology and velocity extent of the structures in 
the ISM. In the analysis described above, we have shown 
that the measured linewidths are rather insensitive to the 
location in PPV space where they are measured. We now 
consider whether the precise shape of the defined structures 
affect the linewidth measurement by analysing regions with 
arbitrary shapes. We consider two simple cases, one in which 
the linewidths are measured in spherical regions, and one 
using cubic regions from the PPV cube. 

Figure [3 shows the result of this exercise, comparing the 
spherical and cubic derived linewidths to a representation of 
the original dendrogram derived relationship from the N2H + 
observations. To portray the original measurement, we em- 
ploy a "posterior prediction" from the Bayesian regression 
of the N2H + dataset shown in Figure [H Posterior predic- 
tion is a straightforward way of comparing the Bayesian fit 
with new data. In other words, it can "predict" new (or fu- 
ture) measurements b ased on the model derived from th e 
existing data (see e.g. I German et all 2004; Kruschke 2 01lh . 
From the Bayesian posterior, at each MCMC step we can 
draw a large sample of a for chosen values of R. Due to the 
uncertainties and the dispersion about the regression line, 
the posterior defines a range in possible a values. The 95% 
HDI of the predicted a values from the dendrogram derived 
a — R relationship is shown by the green lines in Figure 



4 In the limit of completely random and incoherent motions, 
there would be no correlation between linewidth and size, or 
b w 0, as discussed bv lShettv etafl (1201 ll . see their Fig. 13). 



and the short horizontal lines show the mean value. Subse- 
quent data that fall within the extent of the vertical lines 
are considered to be consistent with the Bayesian inferred 
a — R relationship. 

The blue circles in Figure [71 show the measured 
linewidths in randomly extracted spherical regions in the 
N2H + PPV cube. The sizes of the spheres are determined by 
the chosen value of the radius. There are up to three spheres 
for each size considered, but as the measured linewidths of 
each one are nearly identical, it is difficult to distinguish the 
multiple structures. The red squares are measurements ex- 
tracted from cubic regions. The sizes R of the squares are 
computed in the same manner use d to define the sizes of 
the dendrogram structures, R — yj 'A/ty = L/y/ir, where L 
is the length of one side of a face of the cube lying on the 
2D plane spanning the spatial dimensions of the PPV cube. 
The extent of the spheres and cubes along the velocity di- 
mension of the PPV cube is rather arbitrarily chosen. The 
measured velocity dispersion simply scales with the velocity 
extent over which it is measured. For the cubes considered in 
Figure [7] we simply set the extent in the number of channels 
to be equal to the number of spatial pixels defining the cube. 
Similarly, for the sphere the number of spatial pixels defining 
the "radius" sets the extent in velocity. Had we increased the 
extent in velocity, the cube would instead be a cuboid (a 3D 
rectangle) , and the sphere would be a cylindrical region ori- 
ented parallel to the velocity axis, with semi-spherical ends. 
Yet, the size would remain the same, due to its working def- 
inition. Clearly, for the cubes and spheres in Figure [7] as we 
have defined them the measured linewidths are dependent 
on the spectral resolution of the dataset. 

In general, the precise manner in which the sizes are 
defined will affect a ny inferred a — R relationship, as dis - 
cuss ed previously bv Ballesteros-Paredes &; Mac Lowl ([2002) 
and I Shettv et al . (2010). Therefore, such caveats must not 
be ignored in interpreting any apparent relationship. Bear- 
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ing these caveats in mind, Figure still reveals a few inter- 
esting properties. 

First, the a — R relationship measured using spheres 
is similar to that derived from dendrograms on the smaller 
scales where R ^ 5 pc, but diverges at larger scales. Given 
our definition of a sphere and cube, along with our uncer- 
tainty estimates of the size (2 pc) and dispersion (2 km 
s _1 ), a scales (almost exactly) linearly with R. The a — R 
relationship recovered from the spheres falls within the 95% 
HDI of the Bayesian predicted a values. Given the Mopra 
resolution, it is therefore difficult to distinguish a a — R re- 
lationship obtained using the projected morphology of the 
structures, where the most likely index b < 1, from that 
derived through randomly extracted spherical regions as we 
have defined it here, which results in an index b «1. Indeed, 
that a linear relationship falls within the 95% HDI of all 
the tracers (see Tab. [2]) already alludes to this conclusion. 
Within the uncertainties, the indistinguishability between 
the a — R trends measured using the projected morphology 
of structures or using simple spheres indicates that there is 
some class of shapes - which relates the spatial and spectral 
extent of clouds - which would all recover the underlying 
linewidth-size relationship. This may be why previous ef- 
forts using different structure identificat ion schemes of GC 
clouds have found similar results (e.g. iMivazaki &; Tsuboil 
l200d : lOka et al.lfl998l l200lh . 

Second, the cubic a — R relationship systematically re- 
covers higher linewidths than the spherical case, though the 
trend is still linear. These linewidths occur beyond the 95% 
HDI, suggesting that linewidth measured in cubic region 
cannot reproduce the original dendrogram derived a — R re- 
lationship. The offset by itself is not surprising, since the 
cubic regions contain a larger number of zones at velocities 
towards the edges of the defined region, thereby increasing 
the number of velocities with large differences from the mean 
value, correspondingly increasing the dispersion. But, taken 
together with the offset between the spherical and dendro- 
gram defined regions, there is an unambiguous decrease in 
linewidth from the cubes to the dendrograms. The variation 
may be reflective of the role of turbulence, and the coherent 
velocities it generates, in sculpting the dense structure in 
the ISM. 

The shapes of structures formed in the ISM are likely to 
be governed by the nature of the turbulence. We have shown 
that it would be difficult to distinguish between the a — R re- 
lationship derived considering spherical structures or a more 
standard approach using the morphology of the structures. 
But, as shown in Figure [71 the shape of the structures do 
matter, and arbitrary shapes with similar sizes cannot gen- 
erally reproduce the observed trend. Detailed analysis of the 
shapes of the structures, perhaps including more robust defi- 
nitions of the "size," may shed light on this issue. One possi- 
ble avenue is the use of a "volume" instead of size, for which 
3D simulations of the ISM may aid in understanding how 
the shapes and extents of structures in the PPV cube are 
correlated with the true spatial morphology. 



Galactic Center 





R (pc) 



Figure 8. Comparison of the linewidth-size of structures traced 
by N2H+ (black circles) and HCN (red squares) in the Galac- 
tic Centre with 13 CO features from the Milky Way molecular 
cloud Perseus. Lower dashed line is the best fit relationship from 
ISolomon et al . (1987), a = 0.7R 5 . Upper dashed line is the same 
relationship, but where the coefficient is 3.6. 



the star formation process (e.g. [Krumholz & McKed 



2005 



2004 



IPadoan fc Nordlundl l201ll , see iMac Low h Klessenl 
iMcKee fc Ost riker 2007 and references therein). Given the 
uniform trends shown in Figure [2] along with the results 
from local ISM studies, we are in a position to perform a 
comparison between the two environments. 

As a first comparison, and to validate our measurement 
approach, we compute the size-linewidth relationship from 
the dendrograms of 13 CO observations of the Perseus molec- 
ular cloud. Optically thin 13 CO emission traces both the 
diffuse molecular gas, as well as dense filaments and cores, 
within GMCs. It is therefore analogous to the dense gas 
tracers we consider for the CMZ, within which molecular 
gas is pervasive. Figure [8] shows the N2H + and HCN a — R 
relationship in the CMZ, along with the trend measured 
in Perseus. The (x 2 ) best fit line to the Perseus data is 
a — 0.62 i?°' 54 . This rela tionship is similar to the best fit de- 
rived bv ISolomon et al.l ([1987)0, a = 0.7i? 5 , and is shown 
as the lower dashed line in Figure [8] 

Clearly, the a - R trend in the GC is similar to the local 
ISM, in that a comparable power-law index can account for 
the relationship across a range in spatial scale R G 2 — 40 
pc. Yet, there is a systematic enhancement in dispersions, 
resulting in a coefficient that is about five times larger than 
the local ISM value (see Table [1]). The upper dashed line 
in Figure [8] shows the local a- R relationship, enhanced by 
a factor of five. There is a reasonable agreement between 
this line and the linewidth-size tr end measured in the G C, 
in agreement with the analysis bv lOka et al.l (1998, 200l|) of 
individual and distinct CO clouds. 



4.3 Comparison with Milky Way GMCs 

Since the seminal work by [Larson] (|l98lh , the linewidth-size 
relationship of the local ISM has been extensively measured. 
It has played a central role in theoretical efforts to explain 



5 As indicated bv lHever et al.l (|2Q09h , the coefficient must b e re- 
vised to 0.7 from the value reported bv lSolomon et al 1 (I1987T ) due 
to improved estimates of the solar galactocentric radius. 
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4.4 Interpreting the linewidth-size relationship 

The agreement between the a - R trend within a MC and 
between individual MCs has been int erpreted as evid e nce fo r 
the "uni versality" of turbulen ce by iHever &; Bruntl (|2004r ). 
Namely, IHever fe Bruntl (2004) suggest that turbulence ob- 
served i n MCs originates on the largest scale in the ISM 
(see alsolOs senko pf et al.ll200ll : lOssenkopf fc Mac Lowll2002l : 
iBrunt et" al. 2009). In their description, colliding flows in the 
diffuse ISM lead to the formation of MCs. While the MCs 
form at the interfaces of these colliding flows, they inherit 
the turbulent characteristics of the ambient ISM, leading to 
the observed "universal" a - R relationship; smaller scale 
turbulence injection within MCs, e.g. due to internal super- 
novae or stellar winds, appear not to contribute substantially 
to the overall dispersions. 

As we derive the linewidths in both the densest clouds 
and their surrounding medium, we can assess the relation- 
ship of the clouds and its environment. The relatively uni- 
form linewidth-size trend in the CMZ, extending over an 
order of magnitude in spatial scales, may indeed reflect a 
scenario in which the dense clouds inherit the dynamic prop- 
erties of the surrounding gas. Further, our analysis of the 
linewidths irrespective of the position within the PPV cubes 
(in Section T4.1|) advocates for large scale driving as the pri- 
mary source of turbulence. These results suggest that the 
transfer of turbulent energy across spatial scales in the CMZ 
occurs in a similar manner as in the more diffuse and quies- 
cent ISM in the main disc of the Milky Way. 

Nevertheless, the larger coefficient in the a — R rela- 
tionship attests to some fundamental difference in the CMZ 
and the local ISM. The measured star formation rate in the 
CMZ is ^50 times larger than t he rate in the solar neighbor- 
hood (|Yusef-Zadeh et alJl2009h . Correspondingly, there is a 
higher rate of supernovae occurring in this starbursting envi- 
ronment. The variations in star formation activity are tied to 
the differences between the environmental properties of the 
main disk and CMZ. Namel y, the total ISM dens i ty, vertical 
oscill ation period (see e.g. i Qstriker et al.l |2010| ; iKim et all 
l201lh , and perhaps even shear levels set by the rotational 
velocity at a given Galactocentric radius, result in shorter 
dynamical times in the CMZ, and the associated cloud for- 
mation and destruction timescales. In the dense CMZ, the 
crowded sites of star formation and supernovae result in a 
higher frequency of colliding streams, relative to the main 
disc. These dist ributed supernovae may b e a primary source 
of tu rbulence (|Ostriker fc Shettvl l201ll ; IShettv fc Ostrikerl 
2012), which would be external to any given dense cloud. In 
the main disc, GMCs are relatively isolated from each other, 
so that supernovae are not distributed, resulting in a lower 
frequency of colliding shells. Coupled with the larger densi- 
ties and temperatures, the ambient pressure in the CMZ is 
higher than the rest of the Milky Way. The enhanced tur- 
bulent velocity, and correspondingly the larger coefficient in 
the linewidth-size relationship^], is likely associated to the 
higher degree of star formation activity in the CMZ. 

The relationship between turbulent velocities, cloud 
sizes, masses, and the other physical properties of the ISM 
is germane for developing theoretical descriptions of the 



6 See IHever et al.l (2009) for a discussion about the coefficient, 
and its relation to the dynamical state of the clouds. 



star formation process. Recent theoretical efforts have fo- 
cused on the role of turbulence in regulating star for- 
mation (e.g. iKrumholz et al] 120121 : IQstriker etHI l2010l : 
IQstriker Sz Shettvl 12011 ). Numerical simulations of self- 
regulating starbursting systems where distributed super- 
novae driven feedback balances the vertical weight of the 
disk pro duce turbulent velocities in dense gas of order 5-10 
km s" _1 (jShettv fe Ostri ker 2012), similar to the dispersion 
of the (presumably densest) clouds occurring on the smallest 
scales in Figure [2] In these models, the large scale dynamics 
drives the evolution of gas on smaller scales, including the 
velocity dispersion and ultimately star formation in dense 
gas. Quantifying the dynamical properties of the ISM will 
further constrain the numerical simulations and the affil- 
iated description of star formation in similar starbursting 
environments. 

As more observations resolve the structure and dynam- 
ics of a variety of environments, a direct comparison of char- 
acteristics of the local ISM will undoubtedly unfold. Dust 
and/or spectral line observations from lower density tracers 
will allow for estimates of the masse s of the structures in 
the CMZ (e.g. lLongmore et al.ll2012h . Combined with the 
a — R relationship, the full dynamic state of the clouds and 
ambient ISM should emerge. ALMA will be ideal for car- 
rying out such an analysis, and reveal the similarities and 
differences between starbursting regions like the CMZ and 
the more quiescent environments, ultimately leading to a 
thorough understanding of the ISM. 



5 SUMMARY 

We have measured the linewidths and sizes of structures in 
the observed N 2 H+, HCN, H 13 CN, and HCO+ PPV cubes of 
the CMZ. Iso-intensity contours in the PPV cubes identify 
contiguous structures, from which we measure the linewidths 
and projected sizes. The ensemble of structures clearly por- 
trays a linewidth a which increases with size R, with a ~ 2 
- 30 km s _1 for structures with R £ 2 - 40 pc. 

Employing a Bayesian regression analysis, we showed 
that the four tracers all produce similar linewidth-size 
power-law relationships. Fits from the tracers with the low- 
est noise levels, optically thin N2H + and optically thick 
HCN, produce indices b G 0.3 — 1.1, and coefficients A £ 
0.9 - 4.7 (Section [3} . 

By randomly shifting the iso-intensity contours within 
which a is measured, we demonstrated that the measured 
a — R trend is rather independent of the presence of 
dense structures. Therefore, turbulent velocities are predom- 
inantly driven on the largest scales in the CMZ, in agreement 
with interpretations of the ISM in the Milky Way disk. That 
there is no strong evidence for a flattening of the a — R re- 
lationship in the CMZ suggests that turbulence is driven on 
scales ^ 30 pc (Section I4.ip . We also discussed that the 
morphology of structures and its velocities is (to some ex- 
tent) likely to be determined by the nature of the underlying 
turbulent but coherent velocity field (Section 14. 2p . 

We compared the a — R trend in the CMZ to that de- 
rived from Perseus. A power law linewidth-size relationship 
with index 0.5, which describes the local ISM well, reason- 
ably captures the characteristics derived in the CMZ. This 
may be suggestive of the similarity in which star forming 
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clouds inherit the turbulent properties of the ambient gas 
out of which they form. However, the relationship in the 
CMZ requires a factor of five enhancement in the coefficient 
relative to the local a — R power-law (Section S3]). The star- 
bursting nature of the CMZ, along with its higher densities, 
temperatures, and the associated pressure increase may be 
responsible for the larger magnitudes in turbulent velocities. 
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